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ABSTRACT 

The objective of this paper is to investigate whether a convective dynamo can 
account quantitatively for the observed lower limit of X-ray surface flux in solar- 
type main sequence stars. Our approach is to use 3D numerical simulations of 
a turbulent dynamo driven by convection to characterize the dynamic behavior, 
magnetic field strengths, and filling factors in a non-rotating stratified medium, 
and to predict these magnetic properties at the surface of cool stars. We use 
simple applications of stellar structure theory for the convective envelopes of 
main-sequence stars to scale our simulations to the outer layers of stars in the 
FO-MO spectral range, which allows us to estimate the unsigned magnetic flux 
on the surface of non-rotating reference stars. With these estimates we use the 
recent results of Pevtsov et al. (2003) to predict the level of X-ray emission 
from such a turbulent dynamo, and find that our results compare well with 
observed lower limits of surface X-ray flux. If we scale our predicted X-ray fluxes 
to Mg II fluxes we also find good agreement with the observed lower limit of 
chromospheric emission in K dwarfs. This suggests that dynamo action from a 
convecting, non-rotating plasma is a viable alternative to acoustic heating models 
as an explanation for the basal emission level seen in chromospheric, transition 
region, and coronal diagnostics from late-type stars. 

Subject headings: stars: magnetic fields — activity 



1. Introduction 



Understanding the origin of magnetic activity in stars has been an important research 
area in astronomy and astrophysics for many years. Imaged solar observations show a clear 
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link between magnetic fields and the formation of heated plasma in the chromosphere, tran- 
sition region and corona, especially in active regions (see Fisher et al. 1998). Active regions 
are believed to form from loops of magnetic flux that emerge from the base of the solar con- 
vection zone (see e.g., reviews by Fisher et al. 2000; Fan 2004 ^). The large-scale magnetic 
field on the Sun is believed to originate via a global-scale dynamo, powered mainly by the 
velocity shear (differential rotation) between the convection zone and the radiative interior. 
It is often assumed that the amount of differential rotation increases as the rotation rate 
itself increases. For example, theoretical studies of differential rotation mechanisms such 
as the "A effect" (Riidiger 1989) find that the amount of differential rotation generated by 
Reynolds stresses is proportional to the rotation rate. 

Skumanich (1972) was among the first to propose an observational connection between 
rotation rate and the level of activity in stars through dynamo action. Noyes et al. (1984) 
went on to show a clear correlation between the Ca II H-l-K surface fiux and rotation rate, 
arguing that the correlation is made better with the inclusion of a convective turnover time 
via the Rossby number. Vilhu (1984) finds a similar result for chromosphcric, transition 
region, and coronal emission diagnostics. Basri (1987) shows that the activity of dwarf and 
RS CVn stars both correlate well with rotation rate, but with a slight offset between the two 
classes of stars, while the two classes appear to follow the same relationship when the rotation 
period is divided by the convective turnover time (see also Johns-KruU, Valenti, & Linsky 
2000) . These studies provide strong observational support for the existence of the so-called 
a-D, dynamo (e.g., Durney & Latour 1978), which relies on the interaction of convection 
and rotation — specifically, differential rotation (the "fi" effect) and field regeneration by 
cyclonic motions (the "a" effect) — to generate a magnetic field in a fashion similar to that 
believed responsible for the large-scale solar dynamo. 

For cool stars, a straightforward application of the relationship between stellar rotation 
and magnetic activity described above leads one to conclude that those stars that rotate 
slowly or do not rotate (hereafter collectively termed "non-rotating stars") should show 
little or no activity diagnostics. However, this is not the case. There appears to be a 
lower limit to the emission of activity indicators in the chromosphere, transition region, and 
corona. Given the expected absence of a global dynamo in non-rotating stars, this heating is 
usually attributed to a "basal" non-magnetic mechanism such as heating by acoustic waves. 
Schrijver (1987) first introduced the term "basal fiux density" to refer to that portion of the 
chromospheric and transition region line emission which results from a process other than 
the solar- like cyclic dynamo activity (termed "excess" emission). 
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Evidence for a basal flux and the contention that it is acoustic in origin has been based 
on three observational arguments. First, plots of the Ca II H+K chromospheric flux (as 
observed in the Mt. Wilson HK project, Vaughan, Preston, & Wilson 1978) versus stellar 
color show a lower boundary up against which the stellar observations appear to cluster 
(Rutten 1984; Schrijver 1987; Rutten et al. 1991). Such a lower boundary is also apparent in 
other diagnostics such as Mg II h and k and Si II (Schrijver 1987) as well as in C II, C IV, and 
Si IV (Rutten et al. 1991). A second line of analysis examines flux- flux diagrams for samples 
of stars and seeks to make these as tight as possible by subtracting a color dependent basal 
emission level from one or more of the diagnostics under consideration (Mewe, Schrijver, & 
Zwaan 1981; Schrijver 1983, 1987; Schrijver, Dobson, & Radick 1989; Rutten et al. 1991). 
The basal emission level is suggested to be acoustic in origin, while the "excess" above this 
level is believed to be magnetic in origin. In cases where a basal flux is determined from these 
flux-flux studies, its value is similar to the lower limit of the flux found in flux-color diagrams. 
A third argument for the existence of a basal flux comes from studies connecting rotation 
and activity. Schrijver, Dobson, & Radick (1989) extrapolate the color dependent rotation- 
activity relationship in Ca II H-l-K to zero rotation velocity, flnding a non-zero activity level 
which they identify as the basal flux level. This procedure produces a color dependent basal 
emission similar in strength to that found from flux-color and flux-flux analyses, and this 
basal emission is again argued to be non-magnetic in origin. 

Theoretical estimates of the acoustic heating of chromospheres have been in the litera- 
ture for some time (see review by Schrijver 1995). We mention here only some of the more 
recent results which now show good agreement with the observationally inferred basal emis- 
sion levels. Buchholz, Ulmschneider, & Cuntz (1998) compute time dependent models of 
basal heating for monochromatic acoustic wave models, solving the hydrodynamic equations 
together with the radiative transfer and statistical equilibrium equations in order to predict 
basal emission levels in both the Ca II H+K and Mg II h+k diagnostics. These models 
produce the observed dependence of the lower bound of emission with stellar color, though 
the magnitude of the emission depends rather sensitively on whether complete or partial 
redistribution is assumed for the line transfer. Unfortunately, partial redistribution works 
best for Mg II while complete redistribution gives a closer match to the Ca II observations. 
MuUan & Cheng (1994a,b) extend their acoustic wave models to lower mass column densi- 
ties and find the formation of a relatively cool corona (~ 0.6 MK) when applied to F stars 
and to M dwarfs (T = 0.7-1.0 MK). While this temperature for M dwarfs is lower than the 
coronal temperatures found observationally for both dMe and dM stars (Giampapa et al. 
1996), MuUan & Cheng (1994a) do find that their acoustically heated models can produce 
a surface flux in excess of that observed on the least X-ray active dM stars. Given the 
relative success of wave propagation models at predicting the basal level of chromospheric 
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(and possibly coronal) emission, Judge & Carpenter (1998) examined high spectral resolu- 
tion observations of C II lines in basal level stars and in quiet regions on the Sun to look 
for evidence of upward propagating shock waves as predicted by acoustic heating models. 
No such evidence was found, prompting Judge & Carpenter (1998) to call into question the 
non-magnetic origin for basal emission. 

The role of basal flux in the observed X-ray emission is more uncertain observationally 
compared to chromospheric and transition region diagnostics. The early flux-fiux studies 
(Mewe, Schrijvcr, & Zwaan 1981; Schrijver 1983, 1987; Schrijver, Dobson, & Radick 1989) 
did not attempt to find a basal level in soft X-ray emission. Ruttcn ct al. (1991) did find 
a strongly color dependent lower bound on the X-ray emission of stars in his sample, but 
this was likely due to instrumental detection limits, as the recent volume limited study by 
Schmitt (1997) clearly shows that the lower bound on X-ray emission in cool stars has a very 
weak color dependence at best; which, if it exists, is in the opposite sense to that found by 
Rutten et al. (1991). While their X-ray data probably suffered from detection limitations, it 
should be noted that Rutten et al. (1991) do not find a basal level to the X-ray emission — 
they find that all of the X-ray emission is "excess" in nature, and therefore likely magnetic 
in origin. In a later study, Mullan & Cheng (1994b) argued that the observed X-ray emission 
from late A and early F dwarfs is basal (and therefore acoustic) in nature, while Mullan & 
Cheng (1994a) and Mullan & Fleming (1996) have argued that 90% or more of the X-ray 
emission from dM stars is basal (acoustic) in origin. 

In the more recent, volume-limited study of nearby stars, Schmitt (1997) finds that 
essentially all stars with outer convection zones emit X-rays with characteristics similar to 
the solar coronal emission. It seems, therefore, that all "cool stars" possess coronae. Do 
they then also possess magnetic fields? Schmitt (1997) notes that for the nearby cool stars, 
there is a well defined minimum X-ray surface fiux of log Fx ~ 3.7 (in cgs units) which 
appears largely independent of stellar color. Such a lower bound to the X-ray emission is 
reminiscent of the idea of "basal" surface fiux in chromospheric and transition region emission 
in cool stars, which is usually assumed to be acoustic in origin, as noted above. Recent work 
analyzing both solar and stellar data (Pevtsov et al. 2003) indicates a clear and unambiguous 
relationship between unsigned magnetic fiux and coronal X-ray emission that extends over 
twelve orders of magnitude. In that paper, observations ranging from small patches of the 
quiet Sun to the most active pre-main sequence stars were analyzed. At the small-scale 
end of this study, Kitt Peak magnetograms and Yohkoh SXT images were averaged over a 4 
heliographic degree square area centered at the equator and central meridian of the Sun, with 
reported magnetic fluxes and X-ray fiuences normalized to a single arc-second square area 
on the solar disk, corresponding to the approximate size of a single Kitt Peak magnetogram 
pixel (Pevtsov 2004). This study lends further credence to the idea that at least for X-ray 
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emission, a "basal" heating mechanism may indeed be magnetic in origin, but unhke the 
"excess" emission from global fields, the magnetic fiux has httle or no connection to stellar 
rotation. 

High resolution observations of the Sun indicate a viable magnetic origin for this "basal" 
magnetic activity emission: The existence of a small-scale magnetic component which ap- 
pears to exist independent of the large-scale fields that define the solar cycle, and which 
had been predicted to be generated locally by turbulent convective motions (Meneguzzi & 
Pouquet 1986; Durney, Dc Young, & Roxburgh 1993). Detailed studies of the small-scale 
magnetic field on the Sun from the MDI instrument on SOHO (Title & Schrijvcr 1998; Hage- 
naar, Schrijver, & Title 2003) find strong evidence that small-scale mixed-polarity magnetic 
concentrations observed in the quiet Sun show only a weak solar-cycle dependence. By 
studying the detailed evolution of this field over time scales of roughly one day, the authors 
conclude that this flux is formed by the emergence of small-scale bipoles (e.g., ephemeral 
active regions, Harvey 1993) and dissipated by flux cancellation, where the collision of two 
opposite polarities appears to result in their mutual disappearance (Simon, Title, & Weiss 
2001; Parnell 2002). Title & Schrijver (1998) estimate a fiux replenishment time scale of 
forty hours, while the more recent study of Hagenaar, Schrijver, & Title (2003) estimate a 
replenishment time scale of sixteen hours. In any case, the quiet Sun magnetic fiux — how- 
ever it is generated — appears to correctly predict the X-ray fluence for the quiet Sun in the 
study of Pevtsov et al. (2003). More recent observations (Lin & Rimmele 1999; Khomenko 
et al. 2003; Bueno, Shchukina, & Ramos 2004) using spectropolarimetric instruments in con- 
ditions of extremely good seeing have found evidence of a mixed polarity magnetic fleld on 
even smaller, sub-granulation spatial scales. Taken together, these measurements imply that 
turbulent convection may produce a small-scale mixed polarity fleld over a very wide range 
of spatial scales in the Sun. 

Complementing the observational evidence for small-scale magnetic flelds is a broad base 
of theoretical work. Research in the field of "fast dynamos" (e.g., Vainshtein et al. 1996; 
Tanner & Hughes 2003) indicates that in the presence of wide classes of chaotic fiow fields, 
magnetic energy can grow very efficiently, even in the absence of rotation. In these studies, 
the magnetic induction equation is solved kinematically in the presence of a prescribed fiow 
field. If the fluid streamlines are sufficiently chaotic, the growth of magnetic energy can 
proceed exponentially quickly. Much of the research in this area focuses on highly idealized 
flow flelds (e.g. Vainshtein et al. 1996) and does not account for the role of the Lorentz force 
in limiting the growth of the magnetic fleld. 

The major theoretical breakthrough in this area was made by Cattaneo (1999), who 
demonstrated via a 3D MHD simulation in the Boussinesq approximation that a small seed 
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field embedded in a turbulently convecting, highly conducting plasma can indeed grow ex- 
ponentially until the magnetic energy reaches a level of 10-20% of the kinetic energy of 
the convective motions. Subsequent work by Thelen & Cattaneo (2000) has shown how 
variation of boundary conditions and fluid mechanics parameters affect the development of 
the dynamo, while Emonet & Cattaneo (2001) have focused on computing observational 
diagnostics of the small-scale disordered field. 

Turbulent convective dynamo models of this type must not be confused with global-scale 
dynamo models driven by kinetic helicity in small-scale turbulent motions, which require 
rotation to generate an a effect. An example of such a model is that of Kiikcr & Riidiger 
(1999), who used an o? model to understand magnetic field generation in rotating, fuUy- 
convective naked T-Tauri stars, which are beheved to have little or no differential rotation 
( Johns- KruU 1996). The idea is that hehcal turbulence driven by Coriohs forces can impart 
a systematic twist to field lines encountering the convective eddies, leading eventually to a 
large-scale steady magnetic field. In contrast, the turbulent dynamo that we believe to be 
responsible for the small-scale disordered field in the quiet Sun is unaffected by rotation, 
since the convective turnover time in the topmost layers of the solar convection zone is far 
smaller than the solar rotation period. 

The growing evidence for the importance of small-scale magnetic fields in stellar atmo- 
spheres and the observational evidence for a limit to the X-ray emission in stellar coronae 
leads to the question: Is the turbulent convective dynamo a viable mechanism to generate 
enough X-ray fiux to account for the observational limit? To address this question, in § 2 we 
extend the investigation of field generation by a turbulent dynamo beyond the Boussinesq 
approximation into the "anelastic" regime (Gough 1969; Lantz & Fan 1999; Fan et al. 1999) 
— an approximation to the 3D MHD equations that allows for gravitational stratification 
while still filtering out acoustic waves. Here, we describe the numerical model and present a 
detailed analysis of our dynamo simulations. In § 3, we describe our method of utilizing the 
simulation results to predict the level of X-ray emission in main sequence stars, in § 4 we 
present our results, and finally in § 5 we discuss the implications of our results and summarize 
our findings. 



2. Description of the Model Dynamo 

We solve the 3D system of MHD equations in the anelastic limit (see Lantz & Fan 
1999; Fan et al. 1999 and references therein for a description of the anelastic formalism 
and a discussion of the techniques employed to numerically solve the anelastic system). 
Briefiy, the anelastic approximation results from a scaled variable expansion of the fully- 
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compressible MHD equations. Density fluctuations in the continuity equation are of high 
order in the expansion and are neglected — this effectively filters out rapidly propagating 
acoustic waves, and provides a significant computational savings. The density fluctuations 
are retained, however, in the buoyancy force term of the momentum equation, which allows 
for stratification within the computational domain. Thus, the anelastic formulation is an 
intermediary approximation between a fully-compressible treatment (see e.g., Tobias, et al. 
2001; Dorch 2004) and a Boussinesq model (see e.g., Cattaneo 1999). This treatment is 
well-suited for describing the high-/? plasma of stellar interiors, but cannot directly model 
surface plasmas where P ^ 1 and the acoustic Mach number becomes large. 

In particular, we solve a non-dimensional form of the anelastic MHD equations in a 
horizontally periodic, vertically closed Cartesian domain. Our non-rotating rectangular box 
spans five pressure scale heights vertically, which corresponds to a density difference of 
~ 20 between the upper and lower boundaries. The resolution of the Cartesian domain is 
288 X 288 X 72, which gives an aspect ratio for the box of 4 : 4 : 1. The non-dimensional 
parameters Re and Pr (the Reynolds number and Prandtl number, respectively) arc defined 
as Re = prefVcorwH^ef/ fJ' and Pr = [Ji/ K, and are set to 750 and unity respectively (the Reynolds 
number was chosen to be as large as possible without introducing numerical artifacts). Here, 
i^ref denotes the pressure scale height at the base of the domain, k refers to the coefficient 
of thermal conductivity, and is the coefficient of dynamic viscosity (assumed constant). 
The convective velocity is measured in units of Vconv, where Vconv = ((^refS'-f^ref)^^^, (^ref is the 
non-dimensional super-adiabaticity, and g is the constant vertical gravitational acceleration. 

We begin by dynamically and thermally relaxing a purely hydrodynamic model con- 
vection zone. We initiate convection by introducing a small random entropy perturbation 
within a computational domain with a prescribed background entropy gradient (see Fan, 
Abbett, & Fisher 2003; Abbett et al. 2004), and allow the simulation to progress past the 
time when the convective velocity field attains a statistically steady state. We then set the 
magnetic Reynolds number Rm = VcomHrci/f] of the simulation to 1000, and introduce a 
small, dynamically unimportant magnetic seed field (after t = 164 i/rcf/'^conv of the field-free 
relaxation run). The seed field self-consistently grows and evolves within the computational 
domain as the run progresses (see Figures 1 and 2). The vertical boundary conditions on 
the magnetic field are stress-free and non-penetrating at the bottom, and the field at the 
upper boundary is assumed potential. We note that unlike the simulations of Abbett et 
al. (2004), our potential field upper boundary condition does allow for a small amount of 
horizontally-directed signed magnetic flux to diffuse out of the domain. Also note that our 
"global" diffusion timescale of to = 1000 i^ref/fconv — the characteristic time over which 
a magnetic structure diffuses across H^ef (as specified by our choice of magnetic Reynolds 
number) — greatly exceeds the characteristic convective timescale of the simulation. 
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Figure 1 shows the temporal evolution of the total kinetic (Ek), thermal (Eth), and 
magnetic [Eb) energy fluctuations (each integrated over the entire domain). The quantities 
are normalized by the sum of the three {E^ = E^ + Eth + Eb)- We find that the energy 
of the seed field (initially at 10~^^ of the total kinetic energy) increases exponentially with 
a growth time of ~ 16 i/rcf/'^conv until the magnetic field becomes dynamically important 
{t ~ 350 ifref/^conv), after which the increase becomes approximately linear in time. The 
magnetic energy fully saturates to a time-averaged value of 6.7% of the total kinetic energy 
after t ^ 600 ifref/'J^conv, and ranges from 5.5-8.5% of E^. Figure 2 shows the energies 
on a linear scale — as the magnetic energy increases, the kinetic energy decreases, and 
the thermal energy fluctuations of the plasma increase. Figure 3 is a volume rendering 
of both the entropy perturbations (bottom frame) and magnetic field strength |B| well 
after the field has saturated (t = 710 i/rcf/'^conv)- It is evident that strong magnetic fields 
are concentrated in the narrow, low-entropy vortical downdrafts characteristic of stratified 
convection (particularly in the upper half of the box), and that a greater proportion of the 
total unsigned magnetic fiux resides in the lower half of the box. The net signed fiux in our 
dynamo simulations is always zero, and unless otherwise stated, all magnetic fluxes discussed 
hereafter refer to unsigned magnetic fluxes. We note that even though our simulations exhibit 
the vertical flow asymmetries typical of stratifled convection, and even though we use different 
values of dimensionless fluid parameters and do not achieve the highest numerical resolutions 
explored by Cattaneo (1999), we still find similar and significant levels magnetic energy 
relative to kinetic energy over the simulation domain. Figure 4 shows a "magnetogram" of 
the vertical component of the field near the top of the simulation box (left panel) and near 
the base of the box (right panel). 

To quantify the distribution of unsigned magnetic fiux in the domain (after the dynamo 
has fully saturated), we define a magnetic filling factor 

/o ' /o ' dxdy 

where and Cy refer to the horizontal extent of the Cartesian domain, and the quantity N 
is defined as unity if l-B^I > l-Bzl™* and zero otherwise (li?^!'^"* is a chosen threshold value). 
If we adopt a threshold value of \Bz\^^^ /\BzY^^^ = 0.5, we find that along horizontal slices 
near the top and bottom of the domain — — 2.06 i?ref and zi — 0.45 i^ref respectively — 
the time average of the filling factor < / > varies little, and is quite small: < f{zu) >— 
1.47 X 10~^ and < f{zi) >— 6.09 x 10"^. This indicates that strong flelds are concentrated 
and highly locahzed. Here, jS^I"^^ refers to the maximum value of \Bz\ at Zu sampled 
every 5 Hj-d/vconv in the time interval between t = 600-850 i^rcf/'^conv If we choose a cutoff 
of \Bz\'^"-'^ /\Bz\'^'^^ — 0.05 we find a much larger disparity for the same horizontal planes: 
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< f{zu) >— 0.00473 near the surface, and < f{zi) >— 0.149 near the bottom. This suggests 
that the weaker field is more evenly distributed, particularly in the lower half of the domain. 
Figure 5 shows the time averaged filling factor as a function of \Bz\™^ /\Bz\"^^^ along both 
slices. The dashed and dotted lines denote the horizontal planes near the surface and closer 
to the lower boundary respectively (these planes are the same as those shown in Figure 4). 
The solid line denotes the time averaged volumetric filling factor (the volumetric analogue 
to the area filling factor with 15^1™*^ defined as the maximum value of \B\ in the volume). 

A topic of great current interest is the dependence of dynamo behavior on the magnetic 
Prandtl number Pm, which is the ratio of kinematic viscosity to magnetic diffusivity, (or 
the ratio of the magnetic Reynolds number to the viscous Reynolds number, R^/ R^,). Con- 
sidering only the processes of Coulomb collisions in a fully ionized upper stellar convection 
zone, one can show Pm ~ 1.3 x 10~^T^/n, where T is the temperature in K, and n is the 
number density of hydrogen atoms in cm~^. Substituting typical values for T and n, one 
finds Pjn « 1. 

While a complete investigation of the dependence of our own convective dynamo on 
is beyond the scope of this paper, as the simulation described above took 5-6 CPU months 
of computing resources, we have performed several exploratory calculations to get a rough 

idea of how the convective dynamo effectiveness depends on Pm- First, we have used the 
initial convective state of the dynamo run of this paper to explore how P„i affects the early 
growth phase of our convective dynamo. The result is that the growth rate decreases rapidly 
for Pm < 2/3, but shows signs of saturation for P.^ ^ 2. Second, we have used the saturated 
dynamo state of this paper as the initial condition for a series of simulations in which the 
same viscous Reynolds number is used, but where the magnetic Reynolds number has been 
changed, resulting in differing values of P^. The convective dynamo then begins to relax 
to a new state that refiects the changed value of P^. It appears that both the magnetic 
energy and the unsigned magnetic fiux drop significantly for P^ < 2/3, but shows signs of 
saturation for P^ > 2. Beyond this, for a fixed value of Pm, we also find evidence for a 
dependence on the magnetic or viscous Reynolds number due to the fact that we cannot 
simulate the full inertial range. We expect, as per the discussion below, that there will be 
an asymptotic limit to the kinetic and magnetic energies as the Reynolds numbers increase. 

From this behavior, we conclude that for values of Pm much smaller than 2/3 (for the 
range of magnetic and viscous Reynolds numbers that we can simulate), the convective dy- 
namo is ineffective at generating magnetic energies that approach those of the kinetic energy 
of convection. We also conclude that values of Pm ~ 1 — 2 can generate significant amounts 
of magnetic energy and flux, with levels that are within a factor of ~ 10 of equipartition 
with kinetic energy. We speculate that larger values of Pm may to some degree increase the 
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effectiveness of the dynamo beyond the cases we can investigate here, but it seems unhkely 
that the dynamo can generate more magnetic energy than the kinetic energy equipartition 
value. Further, the simulations of Longcope, McLeish, & Fisher (2003) , in which a closed 
magnetic ring is stretched by random convective motions, essentially has Pm = oo since the 
magnetic resistivity is zero in calculations like these. Yet these simulations still attain a 
finite amount of magnetic energy. In the large Pm limit, the field strength is most likely de- 
termined by a dynamic force balance between ram pressure gradients and magnetic tension 
forces, rather than by a balance between field amplification via stretching and field decay by 
magnetic diffusion, as occurs in the kinematic limit in many dynamo simulations. 

Given the dependence of the convective dynamo on P^ that we seem to find, the small 
value of Pjn set by coUisional processes in stellar convections zones suggests that convective 
dynamos are extremely inefficient. However, a number of observational and theoretical 
considerations suggest the appropriate value for is in fact much closer to unity in real 
stellar convection zones. We discuss these here. 

First, magnetic elements on the Sun are observed to move, disperse, and cancel in re- 
sponse to granular and supergranular motions (Simon, Title, & Weiss 2001; Parnell 2002), 
indicating that long before the tiny molecular dissipation scales are reached, motions as- 
sociated with convection can effectively disperse a magnetic element on much larger size 
scales. This means that the effective magnetic diffusivity is much larger than the molecular 
value. A similar argument applied to the observed correlation length and correlation time 
of convective motions also leads to a kinematic viscosity much greater than the molecular 
value, and qualitatively close to the effective magnetic diffusivity - implying a value of Pm 
near unity. 

A related theoretical argument for Pm ~ 1 comes from two observations of our own 
simulations: (1) Both the magnetic and viscous Reynolds numbers of the real physical system 
we arc attempting to model are far greater than those in our code. Therefore, at the smallest 
scale we can resolve in our code, magnetic fields and momentum are still well within the 
inertial range of the turbulent motions excited by convection, and therefore both momentum 
and magnetic fields will be diffused by eddies at the resolution scale and smaller. Applying 
a common eddy diffusivity for both quantities leads to a value of Pm ~ 1. A pitfall to this 
argument would occur if one could show that what occurred on sub-resolution scales could 
strongly affect what happens on the largest scales. This leads to our second observation: (2) 
Both the magnetic and kinetic energy spectra exhibit peaks at large scales, with significantly 
reduced energy at the smallest scales in the simulation, especially near the bottom of the 
simulation where the magnetic energy density is highest. Given this, it seems unlikely that 
details of the dissipation physics occurring on unresolved scales will have a strong effect on 
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the macroscopic properties of the dynamo simulation we have computed using = 4/3. It 
is interesting to note that our conclusion that both magnetic and kinetic energies peak at 
large scales differs from the results of dynamo calculations using the "Kazantsev" approach, 
in which either the velocity (in kinematic models, e.g. Boldyrev & Cattaneo 2004) or an 
external force field (in dynamic models, e.g. Schekochihin et al. 2004) are driven by an 
assumed temporal white noise spectrum. Because our driving motions are computed self- 
consistently from the convection itself, and not from ad-hoc assumptions about a forcing 
term, we feel we have applied our results in a physically self-consistent manner by assuming 
a value of Pm near unity. 

There are other theoretical arguments against using the molecular diffusivities. Long- 
cope, McLeish, & Fisher (2003) argue that the presence of fibril magnetic fields in the 
convection zone plasma results in an effective viscosity that is at least as great as the turbu- 
lent viscosity, if not larger. Abbett et al. (2004) show that magnetic fiux tubes are dispersed 
in a model convection zone with an effective magnetic diffusivity given by a turbulent eddy 
diffusivity. 

Taken together, all of these results argue for a common eddy diffusivity for both processes 
and therefore a magnetic Prandtl number near unity. We adopt Pm = 4/3 instead of unity 
simply because we want to choose the largest possible values of the viscous and magnetic 
Reynolds numbers that do not result in numerical artifacts. Our preliminary calculations 
indicate no significant differences between using results from P^ = 4/3 versus assuming a 
value of unity. 



Here, we relate the anelastic simulation results described in §2, which are performed 
in dimcnsionless units, to physical units so that the resulting magnetic surface flux can 
be compared with stellar observations. We must therefore assign cgs units to the non- 
dimensional simulation results and connect the velocity and magnetic flclds of the simulation 
to the luminosity, surface gravity, radius, effective temperature and surface density of main 
sequence reference stars. 

The variation of the background temperature Tq{z) and density po{z) with height z 
above the base of our anelastic simulation box is taken to be a polytrope with index m = 
1/(7 — 1) = 1.5 (here 7 refers to the adiabatic index for an ideal gas) and is given explicitly 



3. Connecting the Dynamo Model to the Stellar Envelope 



by 




(2) 
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and 

P0(^)=Pref(^)". (3) 

Here, the subscript "ref" denotes chosen reference values at the base of the box (see Fan, 
Abbett, & Fisher 2003). Values of the thermodynamic variables at the top (or surface) 
of the box are denoted with a subscript "surf". The ratio Tref/Tgurf is one of the fixed 
parameters in the simulation described in this paper, and is set to 7.8, and results in a 
stratified background atmosphere with a base-to-surface density ratio Pref/Psurf of 21.8 and a 
corresponding pressure ratio i-*ref/-Psurf of 170.6. The reference temperature Tref at the base 
is determined from the effective surface temperature T^s by setting Tsm-f = T^s, and then 
using the fixed ratio Tj-cf/Tsuri to determine the actual value of T^cf in degrees K. Similarly, 
the value of surface density from stellar atmosphere models, along with equation 3 is used 
to determine the reference density Pref- 

In the anelastic simulation, the unit of length is the pressure scale height at the base 
of the domain, Hj-^f = RTj-d / (flgrcf) , where R is the universal gas constant, /i is the mean 
molecular weight, and gj-d is the value of the gravitational acceleration at the base of the 
simulation box. For simplicity, we assume that jl = 0.5, which corresponds to a fully 
ionized hydrogen gas. Near the surface of the later type stars we consider, the gas will be 
nearly neutral, however. Our neglect of the effects of variable ionization in the background 
atmosphere is a limitation of our current calculations that could be improved in future work. 
We assume that the values of g^ei and gsurf, the surface gravitational acceleration, are equal, 
since the depth of the box is quite small compared to a stellar radius for the cases we consider. 
The depth of the box is given by 

i:. = (m+l)^H-^i/ref, (4) 
-tref 

and the horizontal dimensions of the box, a-nd Cy are four times larger than C^- 

The unit of velocity in the simulations is defined in terms of a dimensionless entropy 
gradient (see § 2). To convert the velocity units to physical values, we note that in our 
simulation, the dominant form of energy transport from the base to the top is thermal 
convection. We therefore require that the vertical energy flux match that needed to carry 
the stellar luminosity in a particular stellar model. In other words, for a given star with 
luminosity and radius R^, the energy flux carried by convection in the outer convective 
envelope must be given by 

We have assumed that throughout the convective envelope, the fraction of energy carried by 
radiative diffusion is far less than that carried by convection, which is a reasonable assump- 
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tion for the stars we consider here. This assumption will break down near the photosphere, 
where radiation begins to dominate. 

To estimate -Pconv in a consistent way that apphes to each of the different stellar types we 
consider, we use a simple form of mixing length theory to obtain an approximate reference 
value of the convective velocity at the depth corresponding to the base of our simulation 
domain. Here, we adopt the mixing length formulation of Mihalas (1978). Combining his 
equations 7-68 and 7-69, and explicitly evaluating Cp as 5R/{2fi), we find 



^conv = — P^^conv > (6) 

a 

where p is the background density of the convective envelope, Vconv is the characteristic 
convective velocity, and a — I /Hp is the ratio of the "mixing length" to the pressure scale 
height. While mixing length theory is a crude approximation, it has been shown to provide 
roughly the right velocity amplitude in simulations of realistic surface convection (Abbett 
et al. 1997), who also find best results when a ~ 1.5. Our use of equation 6 is to take the 
stellar value of -Fconv and the background stratification of p to estimate the value of fconv 
at the base of our simulation (which we will take as the reference value t'ref), and to then 
assume that our simulation velocity is measured in units of this value. If desired, one can 
work backwards to determine the corresponding entropy gradient, but this is not necessary. 
We have used equations 5 and 6 to estimate the convective velocity in the layers just below 
the solar photosphere and find values of ~ 3 km s~^, in approximate agreement with much 
more detailed models of the solar interior (see e.g., Abbett ct al. 1997; Stein & Nordlund 
1998; Asplund et al. 2000; Samadi et al. 2003), and with observations of convective velocities 
from granulation (Hirzberger et al. 2001; Roudier et al. 2003). 

Because the depth of our simulation box is small compared to the stellar radius in all the 
cases we consider, we are justified in ignoring the change in radius with depth and assume 
that Fconv is uniform with depth within the box. This means that f conv at the base of the 
box (i.e., Vref) can be found from equations 2, 3, 5, and 6, and is given by 

( m 

where 

''■"'•=(io^)'- 

Note that equation 7 is extremely insensitive to the precise formulation used in the mixing 
length theory. For example, if the coefficient "10" in equation 6 was changed by a factor of 
ten, it would result in only a factor of 2.15 difference in the resulting value of v^ef- 
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The magnetic field strengths of the simulation are given in units of 

By-ei = (47rpref)^t'ref • 



(9) 



To convert to physical units (G), we simply substitute the expression for v^ef (equation 7) 
into equation 9: 



5, 



ref 



712 



5RI 



r)6 

r surf 



Tr 



ref 



surf 



(10) 



The scaling relationships given here allow us to estimate the amount of magnetic flux and 
the mean fleld strength near the top of our anelastic simulation of the turbulent dynamo. We 
note that although the anelastic formulation is not well suited to modeling the surface layers 
of stellar atmospheres, and that the convection in our simulations is not driven explicitly 
by radiative cooling in the surface layers, our simulation remains appropriate for our study 
since we are only interested in obtaining an estimate of the total magnetic flux threading the 
top of the convective envelope — not the detailed dynamics or distribution of magnetic flelds 
across a stellar photosphere. Table 1 lists the parameters used to calculate the scaling of the 
simulation results for a sample of main sequence stars ranging in spectral type from FO to 
MO. All stellar parameters except surface density arc taken from Gray (1992) for spectral 
types F0-K5, and those for MO arc taken from Rcid & Hawlcy (2000). The surface densities 
(column 6 of Table 1) are evaluated at the depth where the local temperature is equal to the 
effective temperature (column 2) using model atmospheres from Kurucz (1993). 



4. Results 

We have generated signiflcant magnetic flelds (via a convective dynamo) in a stratifled 
turbulent model convection zone by imposing a dynamically insigniflcant seed fleld on a 
statistically relaxed convective state. Our treatment differs from that of Cattaneo (1999), 
since our domain is highly stratified; however, the energetics of our simulation are similar 
to his results — we find that the magnetic field fully saturates at roughly seven percent 
of the total kinetic energy of the computational domain. We find that the magnetic filling 
factor is small near the surface, and larger deeper in the convective envelope, and that the 
total amount of unsigned magnetic flux is concentrated in the lower half of the simulation 
domain. Strong flelds are concentrated and highly localized, while weaker flelds are more 
evenly distributed. We also note that the correlation time of magnetic flux structures near 
the top of the simulation box is 2-3 ifref/^conv, which corresponds to 1-2 hours for a solar-type 
star. 
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The results of applying the a = 1.5 mixing length scaling to the simulation data are 
presented in columns 2-4 of Table 2: column 2 lists the pressure scale height at the base of 
the domain, column 3 lists the convective velocities at the surface, and column 4 lists the 
total magnetic fluxes of the reference stars, The magnetic fluxes were obtained by first 
averaging the total unsigned flux through a horizontal layer near the top of the simulation 
box {z = Zy) over an interval 250 ifref/i^conv centered at t = 725 ifref/i^conv (see Figure 1). The 
simulated result was then scaled from the horizontal area of the computational domain to 
the surface area of the reference stars. Here, we assume that all areas on the stellar surface 
are uncorrelated and generate magnetic field in the same way and with the same efficiency. 

We estimate the X-ray luminosities of the reference stars through an empirical relation- 
ship between the X-ray luminosity and unsigned magnetic flux (Pevtsov et al. 2003). A fit 
to the data presented in Figure 1 of that paper leads to the relation 

Lx = 0.8940$1-^^^^ . (11) 

The X-ray luminosity, Lx, and surface X-ray flux, Fx, for our sample of reference stars are 
listed columns 4 and 5 of Table 2, respectively. 

We compare the results of our simulation to the volume limited sample of cool stars 
discussed in Schmitt (1997). The stellar data itself comes from both Schmitt, Fleming, & 
Giampapa (1995) and Schmitt (1997). These papers provide the X-ray luminosity and cither 
the absolute visual magnitude. My. or the spectral type and My for each star. The observed 
B — V Ys, also given for about half the stars in the sample. The spectral type (if given) or 
My is converted to a S — y index (if needed) and a stellar radius using the main sequence 
cahbration given by Gray (1992) for stars K7 and earlier, or that given by Reid & Hawley 
(2000) for stars MO and later. In the case of the M dwarfs, Reid & Hawley (2000) do not 
give a direct spectral type — My (or B — V^ calibration, so we used the data provided 
in their appendices to construct one. We binned the stellar data for the 8 pc sample into 
spectral type bins (MO, MO. 5, Ml, Ml. 5, . . . ) and averaged the values of My and B — V 
in each bin. These averages were then fit with a 5th-order polynomial to construct our final 
M dwarf calibration of My and i? — V" as a function of spectral type. Using the published 
X-ray luminosities and our calculated stellar radii, we then found the X-ray surface flux for 
each observed star (shown as a filled circle in Figure 6) . Our predicted model surface fiuxes 
calculated using a mixing length parameter o; = 1.5 are shown with asterisks connected by 
a dashed fine in the same figure. The gray region shows the range of the predicted surface 
X-ray fluxes if the convective flux calculated by the models varies by a factor of eight (or 
equivalently if the velocity scaling varies by a factor of two) above or below that calculated 
using a = 1.5. 

It is reasonable to assume that the heating mechanism that produces X-ray emission 
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also produces emission in chromospheric and transition region diagnostics, contributing to 
the basal emission level observed in these lines. Without a detailed magnetic heating model, 
we can only explore this issue using fiux-fiux relationships between X-ray emission and 
diagnostics such as Mg II emission. Many investigators have found good correlation between 
these diagnostics, and here we use the flux-flux relationship 

based on Tabic 4 of Rutten et al. (1991) to determine an estimate of Mg II emission for our 
reference stars. 

Rutten et al. (1991) do not flnd a basal emission level in their X-ray observations, but do 
find one in the Mg II emission which they attribute to acoustic heating. Subtracting off this 
basal emission from observed Mg II emission for all stars then improves the correlation of 
Mg II with X-ray emission. The assumption is that this excess Mg II emission is all magnetic 
in origin, as is the X-ray emission. Since the emission we calculate is explicitly magnetic in 
origin, we use the relationship between excess Mg II emission and X-ray emission (column b 
of Table 4 in Rutten et al. 1991). We note though, that using the total flux- flux relationship 
given in Rutten et al. (1991) increases the predicted Mg II emission by only ~ 0.2 dex, which 
is well within the range of variation shown in the plots for different choices of the mixing 
length in our magnctoconvcction simulations. Since our predicted X-ray surface ffuxcs do not 
depend strongly on stellar B — V, it is not surprising that we find a similarly weak dependence 
on stellar spectral type of the resulting Mg II surface flux: we flnd log(F]y[g jj) = 5.50 at FO 
and 5.34 at MO. 

At flrst glance, these values appear inconsistent with the basal flux observations of 
Rutten et al. (1991) (their flgure lb), where the basal level appears to be log(F|yjg jj) ~ 4.5 
at B — V — 1.5. However, a closer inspection of this flgure (indeed, of all the flux-color 
flgures in Rutten et al. (1991)) shows that the lower hmit to the Mg II emission is deflned 
by giants and subgiants at later spectral types. The dwarf stars are all signiflcantly stronger 
than the apparent lower limit at these later spectral types. In addition, there is no apparent 
piling up of the dwarf star observations against a lower bound as would be expected for 
basal emission (see Schrijver 1995), probably due in part to the small number of dwarf stars 
observed by Rutten et al. (1991) (there are only 13 stars later than B — V = 1.0 in their 
Figure lb). This lack of piling up of the dwarfs is apparent across all B — V (0.3 to 1.5) 
plotted in Rutten et al. (1991). A more complete survey of Mg II emission in K and M dwarf 
stars is provided by Mathioudakis & Doyle (1992). In Figure 7 we plot the Mg II surface 
flux on the dK and dM stars from Tables 1 and 2 of Mathioudakis & Doyle (1992) against 
their observed B — V (we do not plot the dKe and dMe stars as these stars are expected to 
— and do — show much stronger emission than that from basal emission stars). Also shown 
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in Figure 7 is our predicted Mg II surface flux based on our dwarf models from equation 12. 
Our predicted emission levels due to a convective dynamo fall nicely at the lower limit of 
the observed emission in the range that the models are valid {B — V < 1.5). The agreement 
seen in Figure 7 is substantially better than that shown in the acoustic heating models 
of Buchholz, Ulmschneider, & Cuntz (1998), whose predicted surface fluxes are too weak 
(log(F|yjg jj) = 4.6 to 4.8 at B — V = 1.4) and show a much stronger color dependence (for 
B — V < 1.4) than seen in the dwarf star observations of Mathioudakis & Doyle (1992) or 
seen in the dwarf star observations of Rutten et al. (1991). On the other hand, at earher 
spectra types {B — V < 0.5) the acoustic heating models better match the lower bound of 
the Mg II observations than the crude estimates presented here. 

As discussed in § 2, the overall flux achieved in our convective dynamo simulations 
does depend on the value of the magnetic Prandtl number P^. However, the scaling of the 
model to main sequence stars of different spectral types does not — the resulting weak color 
dependence of the X-ray and Mg II emission is independent of Pm, whereas acoustic heating 
models show a strong color dependence. The X-ray data (Figure 6) and the more complete 
Mg II data (Figure 7) show a weak color dependence over the spectral type range in our 
models that is consistent with the convective dynamo predictions and not with the acoustic 
model predictions. 

5. Discussion and Conclusions 

We have performed 3D MHD simulations of a turbulent dynamo in a highly stratifled 
Cartesian domain in order to determine the amount and distribution of magnetic flux gener- 
ated in a non-rotating convective envelope. We use the computed surface flux near the upper 
boundary of the domain and an empirical relationship between magnetic flux and X-ray flux 
(Pevtsov et al. 2003) to determine the lower limit of X-ray emission in main sequence stars. 
Figure 6 suggests that our simple analytic scaling treatment successfully reproduces the ob- 
served lower limit of X-ray flux found by Schmitt (1997) in the range from FO-MO. This 
result suggests that the level of heating by magnetic sources in the coronae of these stars is 
sufficient to account for most if not all of the X-ray flux. 

Stars earlier than spectral type FO are expected to have the thicknesses of their outer 
convective shells rapidly decrease to zero as the surface temperature increases. The low 
surface X-ray fluxes for the three observed stars for which B ~ V < 0.1 are consistent with 
less heating due to a diminished level of magnetic fleld generated by the convective dynamo; 
however, the same argument holds for less heating due to a diminished convective acoustic 
flux. We do not extend our treatment to spectral types later than MO, since we expect that 
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our assumptions of an ideal equation of state and high electrical conductivity will break 
down for cool M dwarf atmospheres. Our treatment implicitly assumes that the Pevtsov et 
al. (2003) relation applies to magnetic flux generated by the turbulent convective dynamo. 
Currently, the best hope for verifying this assumption is forward modeling of coronal heating 
of ensembles of magnetic loops (e.g., Lundquist et al. 2004; Schrijver et al. 2004), as well 
as in understanding the role of sub-pixel magnetic structures on measured magnetic flux 
densities. 

Wc use the flux-flux relation of Rutten et al. (1991) to determine the expected level of 
chromospheric emission in Mg II. For K dwarfs wc flnd good agreement with the observed 
lower limit of Mg II surface flux in the K dwarf sample of Mathioudakis & Doyle (1992), 
suggesting that this observed lower limit is entirely the result of magnetic heating. On the 
other hand, our models predict a very shallow dependence of both X-ray and Mg II emission 
with spectral type, so that they underestimate the Mg II emission at earlier spectral types 
(though our models are a good match to the X-ray emission at these spectral types). This 
suggests some other agent (acoustic heating) is required to produce the minimum observed 
Mg II emission at earlier spectral types, and that this additional mechanism produces rela- 
tively little X-ray emission. At the very latest spectral types, the sharp drop in the Mg II 
surface flux seen in Figure 7 at B — V 1.5 suggests something very different may be 
happening here. We do not expect our models to be valid at these cool temperatures due 
to the relatively low ionization that will result. In addition, the X-ray emission from stars 
observed in Figure 6 does not show this same drop, indicating that the flux-flux calibration 
used to predict the Mg II emission based on the X-ray emission will not hold. We cannot ad- 
dress the level of acoustic heating with our dynamo model. Obviously, if there is convection 
there will be some level of acoustic activity. However, the results we present here provide 
a consistent and viable alternative to acoustic heating for K type dwarfs in the absence of 
a large scale dynamo. At earlier spectral types, our results suggest turbulent dynamos can 
also fully account for the lower limit to the observed X-ray emission, though there may still 
be a substantial acoustic contribution to the coronal emission in earlier type stars. 
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Fig. 1. — Shown are the kinetic (sohd hne), thermal (dotted hne), and magnetic (dashed 
hne) energy fluctuations (normalized by the sum of the three) as a function of time (in units 

of i^ref/t'conv)- 
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Fig. 2. — Same as Figure 1 except on a linear scale. To facilitate comparison, the magnetic 
and thermal energies are vertically shifted. 
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Fig. 3. — Volume rendering of the magnetic field strength |B 
spending entropy perturbations (lower frame) at t = 710 i^ref 



(upper frame) and the corre- 
hcorw (well after saturation). 
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Fig. 4. — The vertical component of the magnetic field at t = 710 -ffref /"^conv along a horizontal 
slice taken at Zu = 2.06 Hj-d near the top of the simulation domain (left frame), and at and 
Zi = 0.45ifref near the bottom (right frame). 
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Fig. 5. — The time-averaged filling factor < / > as a function of for two 

horizontal slices through the domain: near the top (dashed line), and near the bottom 
(dotted line) . Also shown is the volumetric fiUing factor over the entire domain (sohd line) . 
The thin straight lines represent the one count per area and one count per volume baselines. 
See text for details. 
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Fig. 6. — The filled circles represent the X-ray surface flux for each observed star as a function 
oi B — V color. The asterisks represent our theoretical prediction of the lower bound of the 
X-ray surface flux for a choice of mixing length parameter a — 1.5. The gray shaded area 
indicates the amount that these levels can change if the assumed surface velocities change 
by a factor of two. See text for details. 
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Fig. 7. — The filled circles represent the Mg II surface fiux for each observed star as a function 
oi B — V color. The asterisks represent our theoretical prediction of the lower bound of the 
Mg II surface fiux for a choice of mixing length parameter a — 1.5. The gray shaded area 
indicates the amount that these levels can change if the assumed surface velocities change 
by a factor of two. See text for details. 
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Table 1. Stellar model parameters 



Spectral Type 


Te« [K] 


L* [erg s ^] 


R* [cm] 


logg [cm s 2] 


Psurf [g cm 3] 


FO 


6949 


1.84 X 10-^* 


1.05 X 10" 


4.27 


6.71 X 10-* 


F5 


6445 


1.04 X 10^-* 


9.19 X lO^" 


4.32 


1.18 X lO"'' 


GO 


5948 


5.04 X 10^3 


7.52 X 10^° 


4.42 


1.97 X 10"'^ 


G5 


5678 


3.31 X 10^3 


6.68 X 10^° 


4.46 


2.42 X 10"''' 


KO 


5273 


1.75 X 10^3 


5.64 X 10i° 


4.53 


3.02 X 10"''' 


K5 


4557 


6.88 X 10^2 


4.73 X 10^° 


4.60 


5.33 X 10-'' 


MO 


3800 


2.77 X 10^2 


4.32 X 10^° 


4.65 


9.60 X lO"'' 



Table 2. Scaled results for a = 1.5 



Spectral Type 


Hrei [Mm] 


v^ui! [km s 1] 


[Mx] 


Lx [erg s ^] log Fx [erg s ^ cm 2] 


FO 


4.86 


6.67 


6.01 X 10^:^ 


1.86 X 1027 


4.13 


F5 


4.02 


4.99 


4.57 X lO^s 


1.36 X 102'' 


4.11 


GO 


2.95 


3.78 


2.99 X 10^3 


8.34 X 1026 


4.07 


G5 


2.56 


3.32 


2.30 X 10^3 


6.16 X 1026 


4.04 


KO 


2.03 


2.79 


1.54 X 10^3 


3.88 X 1026 


3.99 


K5 


1.49 


1.90 


9.81 X 10^2 


2.31 X 1026 


3.92 


MO 


1.11 


1.23 


7.08 X 10^2 


1.59 X 1026 


3.83 



